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Abstract 



Recently, high resolution angle-resolved photoemission spectroscopy has 
been used to determine the detailed momentum dependence of the supercon- 
ducting gap in the high temperature superconductor Bi-2212. In this paper, 
we first describe tight binding fits to the normal state dispersion and super- 
lattice modulation effects. We then discuss various theoretical models in light 
of the gap measurements. We find that the simplest model which fits the 
data is the anisotropic s-wave gap cos{kx) cos{ky) , which within a one-band 
BCS framework suggests the importance of next near neighbor Cu-Cu in- 
teractions. Various alternative interpretations of the observed gap are also 
discussed, along with the implications for microscopic theories of high tem- 
perature superconductors. 

PACS numbers: 74.20.De, 74.72.Hs, 79.60.Bm 
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1. Introduction 

One of the key issues today in the field of high temperature superconductivity is the 
symmetry of the superconducting order parameter. Penetration depth measurementsB on 
YBCO indicate that the gap has line nodes in momentum space. This information by itself 
is not sufficient to determine the symmetry of the order parameter. Various Josephson in- 



terference experiment sBd on YBCO are consistent with an order parameter which has dx2-y2 
symmetry; some of these directly probe the sign change in the the order parameter under a 
90 degree rotation.! These results are, however, not consistent with c-axis Josephson tunnel- 
ing dataS and other interference experiments.@ The situation is made even more confusing 
by experiments on other high T^. materials: penetration depth measurement^ and tunneling 
datafl on NCCO and tunneling data on the single layer Hg cupratei indicate an apparently 
isotropic gap. 

The improvement in resolution of angle-resolved photoemission spectroscopy (ARPES) as 
well as the large gaps associated with the high transition temperatures of the cuprates allows 
for the first time the possibility of directly mapping out the momentum dependence of the 
gap. We are indeed fortunate that the (quasi) two-dimensional nature of the cuprates allows 
ARPES to directly measure the spectral function of the electrons. Recently, the Stanford 
groupi presented the first ARPES data indicating significant momentum anisotropy in the 
gap function for the two layer Bi cuprate (Bi-2212), with the gap being large along the CuO 
bond direction (vr, 0) and small, possibly zero, in the diagonal direction (vr, vr) in contrast to 
early ARPES dataO on Bi-2212 which showed no gap anisotropy. More recent work0'0 has 
confirmed the Stanford results and found that the observed gap anisotropy is sensitive to 
sample quality and surface conditionsli2l with some samples showing a much smaller, though 
non-zero,0 gap along the (vr, vr) direction compared with that along (vr, 0). Even in the 
Stanford data the anisotropy decreased significantly with sample aging, probably due to gas 
absorption on the surface.^ It must be mentioned that attempts to see the gap with ARPES 
on other cuprates, particularly YBCO, have not been successful. The reason appears to be 
that in Bi-2212, the cleavage is between the two BiO layers which are van der Waals coupled. 
Thus, the act of cleavage minimally disrupts the sample surface, as opposed to YBCO where 
chain copper — apical oxygen bonds are broken.El 

Recently, our group has measured the detailed momentum dependence of the gap in 
Bi-2212 on very high quality sampleJii employing a spectrometer with improved energy 
resolution (described by a gaussian of standard deviation 8 me V) .111 These results are con- 
sistent with the earlier work discussed above, but with improved resolution, were able to 
demonstrate that the nodes in the gap function were not along the (vr, vr) directions as one 
would have for a gap of dj.2_y2 symmetry, but rather displaced at an angle of 10 degrees to 
both sides of the diagonal direction. The purpose of this paper is to discuss the detailed 
angle dependence of the gap, and see which theoretical models are consistent with such a 
gap function. 

Our main results are summarized below. 

(1) We give a tight binding fit to the normal state dispersion data which reproduces the 
experimentally observed Fermi surface. 

(2) We note that there is no clear evidence for two CuO bands, suggesting that the bilayer 
splitting is either weak or non-existent. 

(3) The occupied area of the Fermi surface corresponds to a hole doping of 17%. 
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(4) We demonstrate that data in the Y quadrant are consistent with the superlattice mod- 
ulation observed in structural studies. 

(5) The observed momentum dependence of the superconducting gap is not compatible with 
pure dx2_y2 pairing. We also argue that the data are not consistent with a dirty d-wave or 
a mixed {s + d, s + id, or d + g) gap. 

(6) The observed gap in the X-quadrant is consistent with an interplay between a pure d^^-yi 
state and the superlattice modulation. However this interpretation is not consistent with 
the data in the Y-quadrant. 

(7) A mixed s ± d gap coming from pairing within a bilayer can fit the data but at the 
expense of having two gaps at each k point. So far there is no evidence in our data for such 
a two-gap spectrum. 

(8) The simplest interpretation of the data is in terms of anisotropic s-wave pairing. An Sxy 
gap function with A(k) = Aq cos{kx) cos{ky) provides an excellent fit to the data. 

(9) pairing would arise within a BCS framework from next-near-neighbor (Cu-Cu) at- 
traction with weak on-site repulsion. 

(10) Several microscopic theories which lead to anisotropic s-wave pairing are discussed in 
relation to the data, including models based on charge transfer, extended saddle points, and 
interlayer tunneling. 

2. Normal State 

We first discuss the normal state data. We will use the notation M = (tt, 0) and Y = 
(tt, tt) with r = (0, 0) where F — M is along the CuO bond direction. The circles and crosses 
in Fig. 1 are the Fermi surface crossings in the ARPES measurements in the normal and 
superconducting states, respectively, and the thick line is the Fermi surface of the tight 
binding fit described below. We used six tight binding functions to fit the normal state 
dispersion data (the latter determined by the peak positions of the ARPES spectra), with 
the functions and their coefficients listed in Table 1. The six conditions used to determine 
the coefficients were (1) two points from the measured energy dispersion in the TY direction, 
(2) two points from the measured dispersion in the FM direction, (3) the measured Fermi 
surface crossing in the MY direction, and (4) the energy at the Y point. The latter is not 
experimentally determined (since it corresponds to an unoccupied state) but is a necessary 
constraint to obtain a good fit. The criteria used was to note that the energy dispersion 
below the Fermi energy looks like that of band theory,^ just reduced by a factor of two, as 
previously noted by Olson's group .0 Therefore, we assumed the same (i.e., half the band 
theory value) to be true above the Fermi energy to obtain the energy at the Y point. The 
resulting energy dispersion fit is shown in Fig. 2. This looks similar to that inferred by 
the Stanford groupHa with two differences. First, the energy of the M point is definitely 
below the Fermi energy (—34 meV), and, second, only one CuO band is seen, not the two 
that might have been expected due to the two CuO planes/unit cell. These differences are 
primarily due to the increased energy resolution of the present experiments. 

This indicates that the bilayer splitting is either very weak or non-existent. Although this 
is consistent with band theory results at the Fermi energy, band theory in addition predicts a 
sizable splitting at the M point related to interactions with BiO bands. Such a splitting and 
the existence of BiO bands below the Fermi energy does not appear to be consistent with our 
data. Absence of bilayer splitting, though, is consistent with theories based on incoherent 
c-axis transport. We contrast this with YBCO, where significant bilayer splitting has been 
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seen in ARPES measurementsE2l (consistent with band theory). This difference is due to the 
presence of chains in YBCO which couple differently to the even and odd combinations of 
the bilayer, which together with the buckling of the CuO layers, causes a sizable splitting.ill 
We also note from Fig. 1 the presence of side sheets in the Y quadrant. The Q vector 
connecting the side sheets to the main sheet is within our resolution equal to the superlattice 
Q vector, 0.21(7r, vr), seen in structural experiments.il The lines going through the side sheets 
are simply displacements of the tight binding fit by this Q vector. We find that all of the 
data is consistent with this picture, and it is not necessary, at least near the Fermi energy, 
to invoke the 2 by 2 modulation proposed by Aebi et alii (we note that the samples of that 
group come from the same source as our own). In principle, more side sheets than two will 
occur, but presumably the intensity of these higher order umklapps are reduced by matrix 
element effects. 

Finally, we note that the area occupied by the observed Fermi surface is equivalent to a 
hole doping level of 17%, the same as that for optimal Tc in LSCO. Using quoted values for 
the stoichiometrics of the cations for our samples0 {Bi2.nSrmCai QiCu2.o50s+x) and the 
above hole count, we estimate an a; = 0.08, which is consistent with a variety of published 
results. 

3. Superconducting State 

The method for extracting the gap has been discussed in detail in our earlier paper .111 The 
Fermi momenta are obtained by finding the minimum separation of the quasiparticle peak 
from the chemical potential (operationally, when the leading edge of the energy distribution 
curve has the maximum slope). Experimentally the peak first approaches the chemical 
potential, then disperses away, with the intensity decreasing rapidly as expected in BCS 
theory. 

The spectral function is assumed to be that of BCS theory with a phenomenological 
linewidth to broaden the delta functions. Fortunately, this linewidth does not enter the fits 
at sufficiently low temperatures since the experimental spectra are resolution limited at 13 K 
where the gaps are determined. The spectral function is multiplied by the Fermi occupation 
factor and then integrated over momenta using the ±1 degree width of the analyzer window 
and convolved with the observed energy resolution function. The momentum integrations 
utilize the dispersion derived from the tight binding fit. The background contribution affects 
the fits shown in Ref. |15| only at binding energies larger than 80 meV and thus is irrelevant 
to the gap determination. 

Besides the overall intensity, representing unknown matrix elements and normalization, 
this leaves only one adjustable parameter in the fit: the absolute value (magnitude) of the 
gap. The resulting gap values are plotted in Fig. 3. The basic point to note is that the data 
are consistent with nodes in the order parameter about 10 degrees away from, and on both 
sides of, the (tt, tt) directions. 

The data differ somewhat in the two quadrants, although it should be noted that the 
data in the two quadrants were taken on different samples. This is because the analyzer does 
not have enough angular range to cover both quadrants in the geometry used. Nonethe- 
less, repeated measurements on different samples yield consistent results for the momentum 
dependence of the gap, with the location of the nodes in the two quadrants being equiv- 
alent. We also note that in any particular quadrant, the measured gap has some sample 
dependence, as found in all ARPES work. 
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4. (i-Wave and Related Models 

In this Section we discuss the extent to which our results are compatible or incompat- 
ible with (i-wave pairing. Some of the experimental evidence for and against c?-wave in 
YBCO was described in the Introduction. From a theoretical point of view it has been 
established^ that if antiferromagnetic spin fluctuations mediate pairing then it must have 
dx^-y^ symmetry. Experimentally the spin fluctuations in Bi-2212 have not been studied as 
intensively as in YBCO, nor are there any phase coherence experiments in Bi-2212. Recent 
NMR experimentsil show a rapidly decreasing Knight shift below T^, giving strong evidence 
for singlet pairing. The low temperature Knight shift and relaxation data suggest a gapless 
state which could arise from d-wave with impurity scattering, but gapless s-wave due to 
magnetic impurities cannot be ruled out. 

Returning to the results of Fig. 3 we see that the data are inconsistent with a simple 
dx-i-y'2 gap, since by definition such a gap must vanish along the (vr, vr) directions where 
course, Bi-2212 has a small orthorhombic distortion, which, even at T^, can lead 
to mixing of d^^-y^ with other symmetries. In YBCO, dx2_y2 can mix with s-wave, but in 
Bi-2212, it can mix only with g-wave (^423 representation of the form (x^ — y'^)xy), and thus 
would still have nodes along the diagonal directions. The difference occurs since in YBCO, 
the orthorhombic axes are along the CuO bond directions, whereas in Bi-2212, they are 
along the diagonal directions (and thus, reflection symmetry about the diagonal directions 
is preserved unlike in YBCO). Of course, s — d mixing can still occur in the non-linear gap 
equations below Tc, but the above symmetry in Bi-2212 in a Ginzburg-Landau approach 
would imply two phase transitions instead of one.@ We also note that an s + d gap simply 
moves the node to one side or the other of the (vr, tt) direction; it cannot lead to nodes on 
both sides of (vr, vr), whereas an s + id gap is nodeless. 

Another alternative to consider is the dirty d-wave state which would produce a gapless 
region around (tt, tt). Even though we do not feel that the error bars on the gap estimation 
allow for this possibility, let us for the sake of argument assume that the small gaps in a 
±10 degree region about the diagonal are all consistent with zero. Within the dirty ci-wave 
theory, when such a large region of gaplessness is produced around (vr, vr) the large gap (near 
the M point here) would also be suppressed. Given that this gap is experimentally found to 
be quite sizable is further evidence against such an interpretation, although we caution that 
the degree of suppression is dependent on whether one assumes Born or unitary scattering. 
A more definitive test of dirty d-wave would be to conduct ARPES experiments as a function 
of disorder and compare to theoretical predictions.^ 
Superlattice Effects 

Next we comment on a possible interplay between a (i-wave order parameter and the 
superlattice modulation in X quadrant. These considerations were motivated by the fact 
that observed nodes in the X quadrant are located at the same points where the main Fermi 
surface sheet and the umklapped side sheets are predicted to cross. It should be noted at the 
outset that to the extent the location of the nodes in the X and Y quadrants are identical 
(within experimental resolution), the superlattice cannot be argued to have anything to do 
with nodes. Nevertheless, if one were to focus only on the X quadrant data, the coincidence 
noted above would require an explanation. 

At this time no firm evidence exists in the data for superlattice effects in the X quadrant, 
primarily because the side sheets are not predicted to be well separated as they are in the Y 
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quadrant. Because of this we do not know yet whether the superlattice modulation which 
exists on the BiO layers influences the CuO electronic structure causing a splitting of the 
sheets at the crossing points. It is possible that the BiO layer primarily acts as a diffraction 
grating for the outgoing photoelectrons (which would be sufficient to explain the effects 
seen in the Y quadrant). Nevertheless, since it will turn out that only a small superlattice 
potential on the CuO planes would give the effect described below, it cannot be apriori 
excluded. 

At the crossing points, then, one has two bands and thus a 4 x 4 secular matrix to 
diagonalize (instead of the simple 2x2 secular matrix in BCS theory). Schematically this 
is given by 



and its eigenvalues determine the quasiparticle energies in the superconducting state. Here e 
and e' are the energies of the two bands measured from the chemical potential, V is the matrix 
element of the superlattice potential mixing the two bands, and A is the order parameter. 
The significance of the parameter s = ±1 will become clear below. It is instructive to 
focus on the analytically solvable case e = e' which corresponds to a locus in k-space where 
contours of constant energy of the two bands intersect. 

First consider a gap function (e.g., an anisotropic s-wave gap) which is even with respect 
to reflections about the (tt, tt) direction. This corresponds to s = +1, since the gap is 
the same for the two sheets at the crossing point. The eigenvalues (for e = e') arc E = 



y {e ± Vy + A'^ . Thus the superlattice potential just acts to "shift" the normal state band 
structure but does not affect the nodes. 

For the case of a dx2_y2 gap, though, the gaps for the two sheets have opposite signs at 
the crossing points, which corresponds to s = — 1 In this case, the quasiparticle energies 
are E = \/e^ + ± V. This is qualitatively different from the previous case since the 
superlattice potential can now affect the locations of the nodes: if V exceeds A, a node in 
E must occur for some value of e where there was none previously. 

To study the second case in more detail we solve the 6x6 secular equation describing 
the main band and the two umklapped side bands using for each band the appropriate tight 
binding energy dispersion in the X quadrant discussed above. In Fig. 4, we show a plot of 
the resulting Fermi surface assuming a superlattice potential of 10 meV (this potential is 
assumed to be the same for all off-diagonal terms; we note that the potential term connecting 
the two side bands to each other is in general different from that connecting the main sheet 
to the side sheets). In Fig. 5, we plot the resulting excitation gap (minimum of E) on the 
middle of the three sheets, assuming a cos^k^) — cos{ky) order parameter. We note the dip 
in the excitation gap near the observed nodes. The finite gap along the V — X direction 
occurs since the middle sheet corresponds to one of the side sheets in this case (the sheet 
nearest the F point along this direction corresponds to the main sheet and has a node; its 
gap is also shown in Fig. 5). Although Fig. 5 has remarkable similarities to Fig. 3, the 
main problem for such a model is that it predicts pure d-wave behavior for the excitation 



( e V A \ 

F e' sA 

sA -e' -V 

VA -V -e j 



(1) 




6 



gap in the Y quadrant with a single node along T — Y, in contrast to experiment. As noted 
above, although the unusual anisotropy of the gap is more obvious in the X quadrant, one 
still sees qualitatively similar behavior in the Y quadrant. 
Pairing within a bilayer 

The final rf-wave related scenario to consider is the interlayer pairing model of Ubbens 
and Lee.ii In this model, it is assumed that the two CuO bands coming from the bilayer are 
degenerate (which is certainly consistent with our data) and that the pairing within a layer 
has dx2-y2 symmetry while that between the two layers has an anisotropic s-wave symmetry. 
Again, one has a 4 x 4 secular matrix to diagonalize for the quasiparticle eigenvalues and 
the result is two excitation gaps of the form d ± s. One of these gaps has a node on one 
side of the diagonal direction, the other on the other side. In Fig. 6, we show a fit of this 
theory to the data in the X quadrant assuming the s-wave component to be isotropic (which 
is probably a reasonable approximation on the Fermi surface) . The smaller of the two gaps 
fits the data quite well, however, there is no clear signature in the data of the larger of the 
two gaps at any point on the Fermi surface. 

More specifically, all the spectra can be fit very well assuming one gap.0 Two gap fits are 
possible, of course, given the finite resolution, but our fits to the spectrum at the maximum 
gap value (where the quasiparticle peaks are the sharpest) indicates that the largest splitting 
between the two gaps that can be accommodated by the data is about 10 meV. This should 
be contrasted with the 18 meV splitting implied by the fit. So we would conclude at this 
stage that the model of Ubbens and Lee is not inconsistent with the data, but that, so far, 
our data show no evidence for a second gap. 
5. Anisotropic s-Wave Models 

We now turn to a discussion of anisotropic s-wave pairing, by which we mean an order 
parameter which has the full symmetry of the lattice. While simple s-wave pairing cor- 
responds to a gap function independent of k, i.e., A(k) = Aq, more general s-wave gap 
functions have the k-dependence of the tight binding functions listed in Table 1. (We will 
not discuss functions, such as | cos{kx) —cos{ky)\, which have a singular k-dependence at the 
node). Examples of anisotropic s-wave pairing are the "extended s-wave" pairing, denoted 
by s* 



corresponding to the first two nontrivial entries in Table I respectively. Note that for the 
purposes of this discussion we shall treat Bi-2212 as being tetragonal ignoring the effects of 
the superlattice. Fig. 3 is clearly consistent with this type of order parameter in that the 
measured gap (a) appears to have an extremum along the diagonal directions which would 
require an order parameter which is even under reflection symmetry about this direction 
and (b) has nodes on both sides of the diagonal directions which is also allowed for this 
symmetry. 

One of the flrst microscopic examples of anisotropic s-wave pairing was the s* solution 
obtained by Littlewood, Varma, and Abrahams,^ in the context of charge fluctuation me- 
diated pairing in a three-band model. Depending on the precise shape of the Fermi surface 




(2) 



and Sxy pairing 




(3) 
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an s* gap has either no nodes (i.e., fully gapped) or two nodes per quadrant.E3 A more gen- 
eral phenomenological analysis was presented by Chen and TremblayS who also discussed, 
among other things, the nodal structure of the s^y state which will be shown to be relevant 
to the ARPES data. 

Before discussing microscopic models, we first describe phenomenological fits to the 
ARPES gap data assuming a sign change in the gap function at the node observed in 
the experiment (so that the gap function is smooth at the node). In Fig. 7, we plot the 
location of the nodes of the s* and s^y gap in the Brillouin zone. We note that the nodes of 
the Sxy gap function are very close to the experimentally observed nodes in the excitation 
spectrum, whereas those of the s* gap function are significantly further from the diagonal 
direction compared with the data. This can be seen more clearly in Fig. 8 where we plot 
these two functions over the observed Fermi surface. In fact, we find, quite surprisingly, that 
a pure s^y gap fits the data very well in the X quadrant as shown in Fig. 9. A least squares 
fit of the form s* + s^y gives only a 2% admixture of s* and one of the form s + Sxy gives 
only a 0.2% admixture of the constant (s) term. Alternatively an s + s* gap can also fit 
the location of the nodes, but at the expense of having too small a gap along the diagonal 
directions (this fit has the two terms comparable in size but with opposite sign). 

Within a simple, single band BCS framework, given an effective electron-electron interac- 
tion — Vk,k'5 the gap function at a temperature T is determined by solving the self-consistent 
equation 

A(k)=i:Vi,.^tanh 

where E{\i.) = y'e^(k) + A'^(k). Here we will "invert" this procedure and infer some features 
of ^,k' which would be consistent with the observed A(k). The simplest potentials which 
give rise to s, s*, and Sxy pairing are attractive on-site, on the near neighbor (NN) site, and 
on the next near neighbor (NNN) site respectively. 

The fact that a constant (s) term is not needed in the fit is significant in that it suggests 
that the large Coulomb repulsion on the Cu site is substantially renormalized in the pairing 
interaction between quasiparticles. We note that although the gap functions discussed above, 
other than the simple s state, satisfy the condition X]k'^(k) = 0, it is not this condition 
which must be satisfied to eliminate the effect of an on-site repulsion in the gap equations. 
Rather, it is the condition X^k A(k)/2i?(k) tanh [E{'k)/2T] = which needs to be satisfied. 
Clearly, a pure s^y function violates this condition which implies such a repulsive term in 
the gap equations must be small.il 

The dominance of the NNN attraction in the singlet channel suggested by the Sxy pairing 
is also very interesting. To avoid any confusion, we emphasize that NNN means the NNN 
Cu site, and not the pairing of holes on neighboring oxygen sites. (The s^y order parameter 
cos{kx) cos{ky) Fourier transformed to real space has peaks in the relative coordinate near 
the (±1,±1) primitive lattice vectors). In most theories proposed to date, NNN pairing 
terms are weak and would play a minor role; in fact, we know of no theory where such terms 
determine the symmetry of the gap. 

We now turn to microscopic models. The charge transfer model of Varma and coworkers 
with s* pairing was already mentioned earlier. The valence fluctuation model of Brandowil 
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also gives rise to anisotropic s-wave pairing where the existence and location of nodes depends 
on the details of the input parameters and on the Fermi surface. 

We finally discuss two other theories: the interlayer tunneling model of Chakravarty et 
al.il and the recent theory of AbrikosovUS based on a weakly screened electron-phonon 
interaction in the presence of extended saddle point singularities. While the physics un- 
derlying the two theories is completely different they share two characteristics: first, the 
dominant part of the interactions are "local in k-space" (for entirely different reasons), and 
second, both theories have the ability of having different order parameter symmetries in 
different materials with comparable unlike in simple BCS models. 

In Abrikosov's theory the poor screening of the electron-phonon interaction causes a 
sharp peak in the forward scattering direction (small momentum transfers). Further the 
local density of states is very large near the extended saddle point singularity or flat band, 
around the M point in Bi-2212, and small away from it. The solution of the gap equation 
thus yields a highly anisotropic gap function which is large near M point and small away from 
it. An additional short range Coulomb repulsion is required to produce nodes as pointed 
out recently by Abrikosov.Ea This adds an overall negative constant to the order parameter, 
thus causing the small order parameter along the diagonal direction to become negative, 
with nodes on two sides of the diagonal in a tetragonal material like Bi-2212. However in an 
orthorhombic material like YBCO which in this theory is modeled as having an extended 
saddle point only along one direction, Abrikosov finds an s + gap with a sign change going 
from the a to the h axis. 

The Chakravarty et al. model is based on ideas of AndersonS that the elementary 
excitations, the holons and spinous, are confined within a plane, and thus the coherent 
transport of single electrons from one plane to another is not possible (this is certainly 
consistent with our observation of only one CuO band). The dominant interaction is the 
tunneling of Cooper pairs between layers which lowers the kinetic energy dramatically and 
thus acts to drive the superconducting transition. This pair tunneling term, which is assumed 
to be completely local in momentum space, leads to the high Tc but does not constrain the 
symmetry of the order parameter. The in-plane BCS interaction, which is a weak additive 
term, is as usual non-local in momentum variables and thus determines the symmetry of the 
gap, with the momentum dependence of this gap being modulated by the k-dependence of 
the tunneling term t^(k). (In Bi-2212 electronic structure calculations suggest that t±(k) 
is large along the CuO bond direction and small along the diagonal direction). In the 
Chakravarty et al. model, the largest of the weak in-plane BCS terms could differ from one 
material to another, thus leading to the possibility mentioned in the beginning of this paper 
of different pairing symmetries for different cuprates. 

6. Conclusions 

In conclusion, we have analyzed various theoretical models in light of our measured gap 
anisotropy in Bi-2212. We find that a simple d-wave model can be made consistent with the 
data in the X quadrant, but not in the Y quadrant, if the superlattice modulation is taken 
into account. We find that the only d-wave model consistent with the observed behavior 
in both quadrants is an interlayer pairing modeS although we remark that so far, we have 
no evidence for the two-gap spectrum required by such a model. The simplest gap function 
which naturally explains the two nodes observed per quadrant is an Sxy order parameter of 
the form Aq cos{kx) cos{ky). Within a one band BCS framework such an order parameter is 
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obtained from a next-near neighbor (Cu-Cu) interaction. Several microscopic theoriesEH'EjE^ 
which give rise to anisotropic s-wave pairing were discussed in this context. 
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TABLES 



TABLE I. Tight binding basis functions used in fitting the experimental energy dispersions as 
described in the text. In this notation, the Y point is (tt, tt). The first column lists the coefficient 
of each term (eV), that is e{k) = ^Cirii{k). 



Cj r]i{h) 

0.1305 1 

-0.5951 5 (cos kx + cos ky) 

0. 1636 cos kx cos ky 

-0.0519 \ (cos 2kx + cos 2ky) 

-0.1117 |(cos 2kx cos ky + cos kx cos 2ky) 

0.0510 cos 2kx cos 2ky 
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FIGURES 



FIG. 1. Fermi surface points measured by ARPES. Open circles were measurements in the 

normal state, crosses in the superconducting state. The thick line represents a tight binding fit 
to the data as described in the text. The thin lines are displaced from this by the superlattice Q 
vector. 

FIG. 2. Energy dispersion from the tight binding fit. The filled circles were data points used 
in the fit (the Y point being obtained as described in the text). Note the extended van Hove 
singularity at the M point. 

FIG. 3. Measured energy gap along the Fermi surface in the X and Y quadrants. The data 
for each quadrant were taken on a different sample. In this notation, degrees corresponds to 
the M — X,Y directions and 45 degrees to the T — X,Y directions, with the angle measured with 
respect to an origin at the X, Y points. 

FIG. 4. Fermi surface determined by the tight binding fit with an additional superlattice po- 
tential of 10 meV. Data points as in Fig. 1. 

FIG. 5. Calculated excitation gap in the X quadrant assuming a d^2_y2 order parameter with 
a coefficient of 12.5 meV and a superlattice potential of 10 meV on the middle of the three Fermi 
surface sheets shown in Fig. 4 (filled circles) and on the inner sheet nearest T (open circles). 

FIG. 6. Fit of the experimental data in the X quadrant to the smaller of the two gaps of the 
Ubbens and Lee model. Filled circles are the data and pluses (+) the two gaps of this model. 

FIG. 7. Plots of the nodes of the s* and Sxy functions in the Brillouin zone (dashed lines). 
The former run diagonally, the latter horizontally and vertically. The solid line marks the Fermi 
surface. 

FIG. 8. Plots of the s* and Sxy anisotropic s-wave functions along the Fermi surface in the X 
quadrant of the zone. 

FIG. 9. Fit of the s^y anisotropic s-wave gap to the data in the X quadrant. Filled circles are 
the data and pluses (-I-) the fit. 
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